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We investigate the sensitivity of a disordered system with diffractive scatterers to a weak external 
perturbation. Specifically, we calculate the fidelity M{t) (also called the Loschmidt echo) charac- 
terizing a return probability after a propagation for a time t followed by a backward propagation 
governed by a slightly perturbed Hamiltonian. For short-range scatterers we perform a diagram- 
matic calculation showing that the fidelity decays first exponentially according to the golden rule, 
and then follows a power law governed by the diffusive dynamics. For long-range disorder (when 
the diffractive scattering is of small-angle character) an intermediate regime emerges where the di- 
agrammatics is not applicable. Using the path integral technique, we derive a kinetic equation and 
show that M{t) decays exponentially with a rate governed by the classical Lyapunov exponent. 
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I. INTRODUCTION 

Quantum manifestations of the classical chaotic dy- 
namics represent a central issue for the field of quan- 
tum chaos. To characterize quantitatively the stability 
of quantum motion, Peres ||] proposed to consider the 
fidelity 



of diffractive nature. For a Gaussian random potential 
assumed here this is equivalent to the condition 



M{t) = |(V'|exp(iH'i)exp(-iffi)|i/')|=^ 



(1) 



where H' differs by a small perturbation from the Hamil- 
tonian H of the system under consideration, and \tp) is 
some original state (wave packet). The quantity (|l|) is the 
probability to return into the state lip) after propagation 
for a time t governed by the Hamiltonian H followed 
by a backward propagation with a slightly perturbed 
Hamiltonian H'. Recently, Jalabert and Pastawski |^ 
argued that for a system whose classical counterpart is 
chaotic the fidelity (|^) will decay exponentially in time, 
with the rate given by the classical Lyapunov expo- 
nent. Their work was motivated by measurements of 
a spin-echo decoherence rate in nuclear magnetic res- 
onance experiments and they gave a name "quan- 
tum Loschmidt echo" to the overlap (^. The paper ^ 
triggered a considerable outbreak of research activity de- 
voted to the sensitivity of quantum chaotic systems to 
external perturbations. In a number of subsequent pub- 
lications |, I 0, |, g, |1^, P p g 0, |1| the 

Loschmidt echo was studied (predominantly by means of 
numerical simulations) for a variety of classically chaotic 
systems and its relation to decoherence problems was dis- 
cussed. These numerical works have confirmed the key 
prediction of |^ that in an appropriate parameter range 
the decay rate of the Loschmidt echo is governed by the 
classical Lyapunov exponent. 

In the present paper, we study the Loschmidt echo (|l|) 
in a different context, namely that of a quantum disor- 
dered system. Specifically, we consider a particle moving 
in a weak quantum random potential. The word "quan- 
tum" means here that the scattering on this disorder is 



(2) 



where Is is the quantum mean free path. This situation 
should be contrasted with the opposite case of a clas- 
sical disorder, for which the disorder-induced contribu- 
tion to the action on a distance ~ d is much larger than 
h, and the representation of propagators in terms of a 
sum over classical orbits (as used e.g. in [||) is justified. 
On the other hand, the standard theoretical tool for the 
quantum-disorder regime is the impurity diagram tech- 
nique. It is therefore natural to attempt to apply the 
diagrammatics to the Loschmidt echo problem. 

We show that indeed the diagrammatic technique can 
be used to calculate the fidelity for short times (where 
it is given simply by the golden rule formula), as well as 
for sufficiently long times (where it decays according to a 
power law reflecting the diffusive character of the classical 
motion) . We demonstrate however that for a sufficiently 
smooth (but still quantum as defined by (||)) disorder 
an intermediate time range emerges, where the diagram- 
matic method breaks down. Using the path-integral ap- 
proach, we calculate the Loschmidt echo in this regime 
and find that it does show the decay governed by the 
classical Lyapunov exponent, which is highly non-trivial 
in view of the diffractive character of disorder. 

The rest of the paper is organized as follows. In Sec. || 
we consider the case of a short-range disorder (d ^ Aq, 
where Aq is the electron wavelength) when the diagram- 
matic calculation works in the whole range of times. 
We identify diagrams corresponding to the short-time 
(golden rule) and the long-time (diffusion-induced power 
law) behavior of the fidelity and evaluate them. Sec- 
tion III, which is the central one for the paper, is devoted 



to the case of a long-range random potential {d ^ Aq) 
when the scattering is of small-angle character. Our con- 
clusions are summarized in Sec. tV. In particular, we 



discuss there a connection between the Loschmidt echo 



2 



problem and a recent activity ||19|, ^ |21|, g3[ devoted 
to quantum interference effects in the regime of quantum 
chaos. 



II. LOSCHMIDT ECHO FOR THE SHORT 
RANGE POTENTIAL 

As discussed in Sec. |, we consider a model of the par- 
ticle moving in a random potential inducing a quantum 
(diffractive) scattering. The Hamiltonians H and H' de- 
scribing the forward and the backward propagation in (|^) 
correspond to two slightly different potentials, 



^2 "2 

jj _ P_ yjj. }J' = — h U' 

2m ' 2m ' 



(3) 



where m is a mass of the particle. In this section we 
study the case of a short-range disorder, with the corre- 
lation length d Ao, which is essentially equivalent to 
a (5-corrclatcd (white-noise) random potential. Thus, we 
have the following expressions for the correlators 

{U{v,)U{v^)) = (C/'(ri)C/'(r2)) - t^6{v, - r^), (4) 

where r is the mean free time (to simplify notations, we 
assume it to be exactly the same for U and ?7'), and v 
is a density of states at the Fermi energy. The difference 
between the potentials 5U = U' ~ U \s characterized by 
another time scale r. 



(5f/(ri)(5[/(r2)) -^<5(ri-r2) 
t:vt 



(5) 
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Clearly, we want to study the effect of a weak perturba- 
tion 5U <^ U or, in the other words, f ^ r. Finally, we 
take the initial state \^) in the form of a Gaussian wave 
packet 



(6) 



where cr 3> Aq = 27r/po is a width of the packet and D is 
a number of space dimensions (we set /i — 1 throughout 
the paper). 

To translate (^ into the diagrammatic language, we 
represent the ensemble-averaged fidelity as an average 
product of four Green's functions 

M{t) = J dri...dr6 (v(ri)G«(ri,r2;t)G-4'(r2,r3;-t) 
X r (r3)V'(r4)G^' (r4 , ; t)G^ (rg , rg ; -t)^ (rg) ) , (7) 



where 



dE 



2tt 



l-iEt 



(8) 



The diagrams are then obtained by connecting four lines 
representing the Green's functions in (|^) via the diffuson 
ladders. At sufficiently short times the leading contribu- 
tion is given by the simplest diagram shown in Fig.^. 

The solid lines in Fig.^ correspond to the impurity- 
averaged Green's functions. 



G«-^(6,p) G^^ ■-(£,? 



R' ,A' I 



£ — ep ± i/2T' 



(9) 



with €p = /2rri^ and the shaded box represents the 
diffuson (see Fig. 



n(Q,^) = 



1 



(10) 



where V = v^t/D is the diffusion coefficient and uq — 
Po/m. Note that this diffuson has a non-zero "mass" 1/f, 
since it represents an averaged product of two Green's 
functions in different potentials, {G^G^ ). We will only 
need this diffuson for zero momentum Q and integrated 
with two Green's functions, therefore it is convenient to 
introduce 



^G«(£-t-^J,p)G^'(e-^J,p)n(0,u;) 



(27r) 



To shorten notations, we will also denote G^'^ (e, p) = 

V'(p)G«-«'(e,p) and G^'-^'(e, p)^ = V*(p)G^'^' (e, p), 
where ?A(p) = (47rcr2)^/4e-(P-Po)''^'/2 is the wave func- 
tion in the momentum representation. 
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FIG. 1: Diagram determining the short-time (golden-rule) 
behavior of the fidelity. 

With the above definitions, the expression correspond- 
ing to the diagram Fig. |l| has the form 



M{t) 



dp de de' 
(2^(2^' 



-i(e-e')* 



X G^(e,p)f(e-e')Gi(e',p) 



(12) 



After a straightforward calculation, we get the following 
result for the fidelity. 



M(t) 



(13) 
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R p + Q/2,e + w/2 



fl(Q,w) = 




j^/ p-Q/2,f-a;/2 

FIG. 2: Dyson equation for the "massive diffuson". The 
dashed hne corresponds to the correlator {UU'}. 



Fig. g do not cancel when all difTuson momenta and 
frequencies are equal to zero. Further, we neglected 
the momentum and frequencies of massless diffusons in 
Eq. (|l5|); the corresponding terms will be smaller by a 
factor f/t^ 1 . 




This is nothing but the golden rule decay induced by the 
perturbation (|5|). 

For long times t ^ f the contribution ( p^ becomes 
exponentially small in view of the massive character of 
the diffusons ( p^ . The long-time behavior of the fidelity 
is however determined by a different diagram shown 
in Fig. ^, with two massive diffusons "colliding" and 
transforming into two conventional, massless diffusons 

(Fig. 3), 



n(Q,L^) 




FIG. 3: Diagram determining the long-time diffusive asymp- 
totics of M{t). 

The diffuson "collision" vertex, conventionally termed 
the Hikami box, is given by a sum of diagrams shown in 
Fig. |, yielding 



x(Q,e-e') = 



1 



Note the unconventional form of the expression for the 
Hikami box: due to {UU') ^ {UU), the diagrams on 



n(Q,w) 




p-l-Q/2,f-|-w/2 



p-Q/2,e-a)/2 



FIG. 4: Dyson equation for the conventional diffuson n(Q, oj). 
The dotted line corresponds to the correlator {UU) = {U'U'). 



FIG. 5: Diffuson "collision" vertex (Hikami box). 

Combining everything, we thus get the following ex- 
pression corresponding to the diagram Fig. ^ 

M{t) = I ^P' f c^e de' dwi duj2 ^-i(^,+^,)t 



i2n) 



3D 



xG>-§,.+ f)G>+§,.-f) 
xn(Q,c^i)x(Q,e - e')|f (e - e')|'n(Q, c^s) 
x<(p'-§,.'-f)G^J'(p' + §,.' + f). (16) 

Performing all the energy integrations, we get the follow- 
ing result: 



M{t) 



xe 



J (27r)3^ 

p'-Pn)V- (W)^ 1 



(17) 



Before writing down the final result, we should be more 
specific about the width a of the original wave packet. 
When it is large compared to the mean free path, a ^ 
vqt, the characteristic deviations |p — Po|, |p' — Pol are of 
the order oil/ a due to Gaussian factors, and we can set 
IpI = Ip'I in the last factor in (p^). In the opposite case 
(7 ^ Vqt the difference |p| — |p'| is of order 1/{vqt) <C I/ct 
and can be neglected in the above Gaussian factors. Thus 
we have 



(15) M{t) = 



2a2 



M(t) 



VQT T(D/2) 
2-KG ipacr)^-'^ yAVt + a 

1 riD/2) f 



2v27r {pQa 



.D-l 



\2Vt 



, a > WOT, (18) 

(T < WOT, (19) 



where T{x) is the Euler gamma function. Remarkably, 
this large-i behavior of the Loschmidt echo is indepen- 
dent of f , i.e. of the perturbation strength. 

Therefore, the long-time asymptotic behavior of the 
fidelity is a power-law decay governed by the diffusion, 
with M{t) proportional to the inverse diffusion volume 
^iff ~ {'Dt)^^'^. Clearly, this result obtained from the 
diffuson diagram technique is not specific for the white- 
noise disorder considered in this section but rather yields 
a generic form of the long-t behavior of M [t] in diffusive 
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t* 

Time, t 

FIG. 6: Schematic representation of the time evolution of 
the Loschmidt echo M{t) for a white noise disorder on a log- 
log plot: A - golden rule exponential decay (^^, B - diffusive 
power law decay (pj|). 



systems. Comparing the results (|l^), ( p^ with contribu- 
tion ( p^ ) , one can determine the time t* of crossover be- 
tween the exponential and the power-law regimes, which 
is larger than f by a logarithmic factor. In particu- 
lar, for a two-dimensional system with a < I we have 
t* = I In £2^. The behavior of M{t) for the case of a 
short-range potential studied in this section is illustrated 
schematically in Fig. |6[ 

It is worth mentioning that our diagrammatic calcula- 
tion bears a certain similarity to earlier studies of inten- 
sity fluctuations and correlations for waves propagating 
in random media ||2j, |2^, |2^. In particular, sensitivity 
of transport quantities to a small change of the impurity 
potential (e.g. due to a displacement of a single scatter- 
ing center) has been investigated |2^. However, the 
Loschmidt echo is essentially different from the quantities 
studied in these papers since it involves propagation in 
two different potentials U, U' already before the ensem- 
ble averaging. Formally, this corresponds to switching 
between the Green's functions G and G' in external ver- 
tices (see Figs. 



of order tt. It is parametrically larger, rtr ~ {d/ X^f'Ts- 
The transport time determines, in particular, the diffu- 
sion coefficient V — v^Tt^/D. 

The condition ^ implies that in the diagrammatic ap- 
proach the leading contribution is given by diagrams with 
non-crossing impurity lines. In particular, Tg is deter- 
mined by the Born-approximation self-energy diagram, 
Ttr is obtained by taking into account the ladder-type 
vertex correction, and the diffusion propagator can be 
calculated by solving the equation corresponding to the 
sum of the ladder diagrams [ p8| . Furthermore, even in 
the ballistic range of frequency and momenta qltt > 1, 
LOTtr > 1, the average product of two Green's functions 
{GrGa) is determined by a sum of ladder diagrams, the 
"ballistic diffuson" . One might thus expect that the di- 
agrammatic calculation of the preceding section can be 
generalized to the case of a long-range disorder. Indeed, 
both the golden rule short-time behavior corresponding 
to Fig. and the long-time diffusion power-law asymp- 
totics determined by the diagram of Fig. || do retain their 
validity for a long-range disorder. However, as we demon- 
strate below, an intermediate time range emerges, where 
the diagrammatic approach is not applicable. The rea- 
son for this is a necessity to average a product of four 
Green's functions describing four electronic trajectories 
propagating close to each other. As was shown in |2^, in 
a certain time range (specified below) these four Green 
functions do not decouple into two (ballistic) diffusons, 
but rather are coupled all together by impurity correla- 
tors into a more complicated object, a "4-diffuson" . In 
view of the failure of the ballistic-diffuson diagrammat- 
ics, we will use the path integral approach developed in 
p3| , p9| . For simplicity, we consider a two-dimensional 
system. 

We begin by defining the disorder correlation functions 
(replacing the white-noise formulas (0) and (||)), 



(C/(r)C/(ri)> 
{5U{v)5U{vi)) 



{U'{r)U'{n)) 
26W{\r- 



M^(|r-ri|), 



(20) 



III. LOSCHMIDT ECHO FOR LONG RANGE 
POTENTIAL 

After having understood the behavior of the fidelity in 
a white-noise disorder, we turn to the case of our main 
interest: a long-range potential with a correlation length 
d ^ Xq. For this kind of potential, the characteristic 
angle of diffraction for each scattering event is small, 
Scj) Xo/d, so that many scattering events are needed to 
change strongly the velocity direction. As a result, the 
motion in such a disorder is characterized by two relax- 
ation times. The first one is the quantum (or, in another 
terminology, single-particle) relaxation time Ts, which is 
the mean time between scattering events. This time de- 
termines the decay rate of the averaged Green's function 
(G«(r,t)) = G^(r,i)e^*/2^=. The second one, the mo- 
mentum (or, transport) relaxation time Ttr sets the time 
scale on which the velocity direction changes by an angle 



Introducing the Feynman path integral representation 
and averaging over the disorder we rewrite the product 
of four Green's functions in (^ as 

(GflXR-i; 1^2, 7")Ga'(R2, R3, —T) 
xGfl'(R4, R5, r)Gyi(R5, R-e, — r)) = 

^ri(T)=R2 /.r3(T)=R3 /.r2(T)=R5 



(21) 



/ri(0)=Ri Jr3(0)=R2 "'r2(0)=R4 
^r4(T)=R6 4 

X / 11-^''^ exp[iS'kin - ^w] 

Jr4(0)=R5 

T 



_m 

>3kin — ~^ 



I dt (f? +rj-r\ 
Jo 



'■4J I 



1 



w 



(^1 



S22 + S; 



33 



tS'4 



+ S12 + S'34 — <5'l3 — — S23 — S24 
— SS12 — SS34 + 6S13 + 3324, 
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5y = / / W{r,{t)-rj{t'))dtdt', 
Jo Jo 

6S^j=[ [ SW{r,{t) -rj{t'j)dtdt', (22) 
Jo Jo 

where the paths ri,r2 correspond to the retarded, and 
r3,r4 to the advanced Green's functions. The path in- 
tegral (|2^) is similar to the one evaluated in Ref. a 
difference being in boundary conditions and in the addi- 
tional terms SSjj in the action induced by the perturba- 
tion SU. As in [E3|, it is useful to perform the change of 
variables, introducing R_|_ — (ri + r2 + rs + r4)/4, R_ — 
(ri -I- r2 - rs - r4), r+ = (ri - r2 -f ra - r4)/2, r_ = 
(ri — r2 — ra -I- r4)/2. The boundary conditions in terms 
of the new variables are as follows. At t = T we have 
R_(T) = and r_(T) = 0, while the integration over 
R+(r) and r+(r) are unrestricted. At t = the integra- 
tion over R±(0) and r±(0) is performed with the weight 



V'(Ri)^*(R3)^(R4)^*(R6) 



1 



X exp 



4R+(0)2 + R^(0)/4 + r+(0)2 + r_(0)2 



2cr2 



+zpoR-(0)}. 



(23) 



The kinetic part of the action reads in the transformed 
variables as 



where 

U ^Uo + SU, 

Uo = 2{F{p+) + F(p_)) - F{p+ + p_) - F{p+ - p_), 
SU = -2SF{p^) + 5F{p+ + p_) + 5F{p+ - p_), (28) 

and we have introduced the functions 



F{y) ^ ^ 


dx 

/ —[W(x,0)^W{x,y)l 
'o ^0 




SF{y) EE 


r ^[SW(x,0)-SW{x,y)], 
Jo ^'0 


(29) 


Giy) ^ ^ 


/■°° dx 5^ 

/ -^W{x,y), 

lo " cly 




SG{y) EE 


Jo V dy-^ 


(30) 



Since the correlation functions W{r), SW{r) decay on 
the scale d, the functions F, G, SF, and SG have the 
following asymptotic behavior: F{y <C d) — — G'(0)y^/2, 
F[y > d) ~ T-\ G(0) = -m^t-Vrtr, and G{y > d) ^ 0, 
and analogously SF(y < d) ~ -SG{0)y'^ /2, SF{y > 
d) ~ T-^, SG{Q) = -m^v^/h,, SG{y > d) ^ 0. Here 



the times Ts and Ttr are defined according to 

1 2 f°° 
— = — / Wir)dr, 

Ts Vo Jo 



(31) 



(*- 



'5'kin = m I dt{ r+r_ -I- R+R- 



(24) 



As was shown in the pairs of variables (R+,R_) 
and (r+,r_) decouple. On ballistic distances <C hi the 
integral over the first pair is essentially of the free-particle 
type, and its saddle-point yields the classical equation of 



motion for the "center of mass" coordinate R 



R+(t)-R+(0) + (R+(T)-R+(0))-. (25) 

After integrating out R+, R_, the action ( p^ is reduced 
to the form 



Skii 



R+(r)-R+(0) ^„ 
m ^ (R-iT) - R-(0)) 



m / dtr+T- 
lo 



(26) 



Since we are interested in the ballistic scales (^ Ztr), it 
is convenient to split R_, r+, r_ into components parallel 
(II) and perpendicular (_L) to the direction of the motion 
R+. Then the disorder- induced part of the action Sw 
depends only on the transverse components R_j_ and 
r±±, which we will denote and p± respectively, 



Sw^ J U{p^{t),p+{t))dt~2 J dtY'^it) 
x{G{p^(t))~5G{p^{t)) + G{p+m, 



(27) 



1 

Ttr 



1 



m-v^ Jo 



dr dW{r) 
r dr 



(32) 



As shown in |29 , these are exactly the expressions for 
the single-particle and the transport times in a long- 
range disorder. The times fj and ftr are defined by the 
equations analogous to (|3l|) , (^2|) but with a substitution 
W SW. 

Taking into account the boundary conditions (|23| ) and 
integrating out the variables R±, we get the following 
expression for the fidelity 



M{T) ^ I ^-P±^e 
' V2^ 



(33) 



The function g entering (|3^) is determined by the po- 
part of the path integral, which can be reduced in the 
standard way to a differential equation: 



d 



dt m dp+dp^ 



+ U{p+,p-) ] g{p+,p-,t) 

Smp^). (34) 



The l.h.s. of this equation reduces to that of Eq.(36) in 
p3[ if the forward and backward evolution are performed 
in the same potential, SU = (the r.h.s. differs from 
p3| because of different boundary conditions) . The pres- 
ence of SU [which enters the "potential energy" U, see 
Eq. (p^)] in (p^) is crucially important: otherwise the 



solution would be simply g{p^, p^,t) = 5{pJ) for any 
t > 0, since the boundary condition in (jsj) is indepen- 
dent on p+ and lA{){p+,()) — 0. After a substitution into 
(^), this would lead to M{t) — 1, which is the correct 
result in the absence of perturbation. 

We have therefore reduced the problem of calculation 
of the Loschmidt echo in the ballistic time range to a 
solution of the kinetic equation (Q) . Let us consider the 
time evolution of the solution of Eq. (jsj). The initial 
value at t — > is g = 5{p-), and as explained above, 
the time evolution is initially determined by the term 51A 
which induces a p+ dependence of the solution. As a 
result, the distribution g becomes quickly suppressed at 
p+ > d because 



SU{p+,0)~^, 



p+ > d. 



Thus for p+'» d Eq. (|3|) reduces to 



d_ 

di 



9 



2 



which gives an exponential decay 
9 - 



Sip-), p+«.d 
6{p-)e-^'/^% p+:$>d 



(35) 



(36) 



(37) 



Therefore, for i 3> the function g remains essentially 
non- vanishing only in the region p-f-,p^ <C d. In this 
region we can expand U up to the leading terms in 
and p_ and rewrite Eq. ( p4[) in the form 



d 



lib T 

+ —P+P- + 



dt mdp^dp^ '''l '''ti 
We have introduced here a new time scale 
3 



■P+ 



3 = 0. 
(38) 



n 



2 

d_ 



dr d /IdWir) 
r dr \r dr 



(39) 



As discussed below, is equal (up to a numerical coef- 
ficient) to the inverse Lyapunov exponent in the corre- 
sponding classical problem, and we will call it the Lya- 
punov time. 

At the early stage of the evolution, characteristic val- 
ues of /9_ are small, and the p\p^^ term in ( ^8|) is small 
compared to the p^_ term. More specifically, at t ~ f ^ we 
have p_ ^ Is/ipod), so that the condition for neglecting 
the quartic term in this time range is II S> h, where 
II = vqTl is the Lyapunov length. We will assume in the 
sequel that this condition is fulfilled ||o). Thus, Eq. (||) 
reduces to 



Performing further a 
i{mvQ)~^d/d4>, d/dp^ 
the following form 



d 



Fourier transformation p+ — s- 
— > imvocf), we cast Eq. (|3J) into 



dp. 



2 92 



(41) 



Remarkably, Eq. (|41|) has a meaning of the Boltzmann 
kinetic equation for the phase-space distribution function 
describing the motion in the transverse direction charac- 
terized by the coordinate with vo(p playing the role 
of the corresponding velocity. This clarify the meaning 
of (f) (and explains the notation): it is the angle the ve- 
locity vector makes with the 1 1 axis. (We remind that 
we are considering ballistic time scales, so that cf) 1) 
The last term in (|4l|) plays the role of a collision integral 
and describes a diffusion process for the velocity angle. 
Solution of this equation is a Gaussian packet 



,P-) 



V37 



exp 



/ 30p_ 



Sp'- 



Transforming back to the variable p^, we get 



2t 



9{P+,P-) = 



where 



1 



exp 



-^-^ + 2^V3^ 



E_(0 
S+(i) 



2vQt / 2t 

\Ttr 

2t 



mvo 



1/2 



1/2 



(42) 



For the phase-space distribution function g{if),p_,) the 
quantities E_ and Y,'^^ play the role of widths of the 
distribution with respect to the coordinate p_ and the 
momentum mvQ(j), respectively. 

To simplify calculations, we will neglect the cross- 
correlations between p+ and p- (this will only influence 
a numerical prefactor in M(t), which is of minor impor- 
tance here) and assume that g has the form 



9iP+,P-) 



1 



■ exp 



p'- 



pI 



(43) 



We will see that this form of g will preserve during the 
further evolution of the distribution. 

The characteristic values of p_ are increasing propor- 
tionally to . Eventually, the neglected third (quartic) 
term in Eq. ( p8[ ) becomes comparable to the fourth one. 
Using the result (^) for the characteristic value 5]_(t) 
of /9_, it is easy to sec that this happens &t t = tl- At 
t > tl the fourth term dominates, and equation (|3J 
takes the form 
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After the Fourier transformation from p+ to 0, the last 
term takes the form —{rn?/Tf^)d'^g/d(tP' and describes a 
diffusion process for the angle 4> with the diffusion co- 
efficient proportional to to the coordinate p-. This is 
exactly the kinetic equation for the disorder-averaged 
distribution function g{(t), p-) of phase-space separations 
between two classical paths ^ . It leads to an expo- 
nential increase of the width of the distribution function 
g{p^,(f>) with a rate given by the Lyapunov exponent 
^ Tj^^ ■ Thus in this Lyapunov regime we have 

S_(i) ~ voTl(^] expic— 



E+(t) ^ Ao ( ^ 



1/2 



exp <; }■ , 



(45) 



where c is a numerical coefficient of order unity. 

Let us emphasize a highly non-trivial character of the 
emergence of a classical kinetic equation. Indeed, we con- 
sider a random potential for which each scattering act 
is of a quantum (diffractive) nature and can not be de- 
scribed classically. It is only after the disorder averaging 
that the classical kinetics is restored. 

The Lyapunov regime breaks down when the width 
E_ reaches the correlation length d, i.e. at r|, = 
{tl/c) ln(ftr/Ttr). This time plays a role analogous to the 
Ehrenfest time but only depends on classical parameters 
(if one assumes that the perturbation 5U is independent 
of h). 

When t > T^, so that p_ :s> d, we can use an- 
other asymptote of the "potential" V({p+, p-)\p_-^d = 
2m^Vn£'^/Tti. This leads to an equation very similar to 
Eq. (EQ) but with ftr replaced by Ttr- Therefore, in anal- 
ogy with ( ^ ) , we have again a power-law dependence of 
the distribution widths, 



S_(<) ~ vot[ — 



1/2 



E+(t) ^ Ao (^) 



(46) 



For times larger than t*^ all the calculations can also be 
performed using the diagrammatic technique for ballistic 
systems. This is because four trajectories, which were 
coupled all together into a "4-diffuson" by disorder cor- 
relations in the Lyapunov regime, split now in two con- 
ventional ballistic-diffusons separated by a distance ^ d 
@,||. 

Having got the solution g of the kinetic equation in 
all the regimes of interest, we can calculate the fidelity 
M{t). Substituting iM) into Eq. (p|), we get 



M(T) 



min(i:+(r),cr) 
max(S_(T), a) 



(47) 



Thus, two additional time scales may become significant 
when we consider the fidelity: Tq-^, which is defined by 







1 2 3 4 5a 5b 5c 6 



TL Tg To-- Toll T-fr 



FIG. 7; Characteristic time scales separating various regimes 
of the behavior of M{t). The regimes 2 (golden rule) and 4 
(Lyapunov) correspond to an exponential decay of M{t), the 
remaining regimes to a power-law decay (see text for details). 



of these time scales with respect to main characteristic 
times (tl, t^,, Tti ) depends on the width a of the initial 
state. As an example, we choose d < a < Itr ■ In this case 
S+ < a for all times t > fs , so that the scale Ta+ does not 
arise, and < Ta- < Th- The order of all characteristic 
scales on the time axis is illustrated in Fig. 0. 

We mention for completeness that there is one more 
characteristic time t„\\ located between t^- and Ttr. At 
this time, the approximation ( p5| ) of straight motion in 
the coordinate R+ loses its validity. This leads to an ad- 
ditional factor a/ maxjcr, in the expression for M{t), 
where 



t 

vt — 

Ttr 



(48) 



characterizes longitudinal fluctuations of R+. There 
is thus an additional crossover inside the ballistic- 
diagrammatics regime, which takes place at a time t^\\ 
satisfying S||(t<^||) = (J- 

We are now prepared to summarize the results of this 
section and to give a list of all the regimes of behavior 
of the fidelity. We have found as much as 6 essentially 
different regimes (as illustrated in Fig. one of them 
(the ballistic diffusion) splits up into three subregimes 
with different power-law behavior. We list the regimes in 
the order they appear as the time increases: 

1. Perfect echo regime, t <C Tg. 

At such short times the perturbation is essentially 
irrelevant, and M{t) ~ 1. 

2. Golden rule regime, t < fsln(cr/d). 

Substituting Eq. (^7|) into (^3|) we get the exponen- 
tial decay. 



Mit) = e-2*/^= 



(49) 



3. Power-law "pre-Lyapunov inflation" regime, 
fs hi{a/d) < t < Tl. 

Substituting Eq. (||) into (E^) we get 



M(t) 



1/2 



(50) 



S+(t, 



(T, and To— defined by I]_(to— ) = cr. Position 



This behavior of the fidelity is related to a power- 
law spreading of classical trajectories in this regime 
due to the "ballistic diffusion" in the perturbation 
potential 5U. 
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4. Lyapunov regime, < t < t^. 

Combining (^5|) and (^) , we get an exponential de- 
cay of the Loschmidt echo determined by the clas- 
sical Lyapunov exponent. 



Ct/TL 



(51) 



5. "Ballistic diffusion" regime, t'^ < t < rtr- 

This regime characterized by a power-law behavior 
of the fidelity is further subdivided into three sub- 
regimes: 

(a) < t < Ta-. Using (^), (|6|) and noticing 
that E_ is still smaller than a in this time 
range, we find 



M{t) 



1 



1/2 



(52) 



(b) Tcr- < t < Ta-\\. Thc Only difference compared 
to the previous case is that now E_ > cr, yield- 



ing 



M(t) 



Pol 



-(-)' 

)Ztr V t J 



(53) 



(c) Tg-ii < t < Ttr- In this regime fluctuations in 
the motion in the parallel direction also be- 
come important, see the text around Eq. (48), 
with the result 



M{t) 



■poll, V t 



(54) 



6. Conventional diffusion regime, t > rti-. 

In this regime the nature of disorder (short-range 
vs. long-range) is irrelevant, and the result ( p^ ) 
derived in Sec. y is applicable, 



M(i) 



_0_7tr 



(55) 



We would like to remind that the ordering of relevant 
time scales in Fig. |^ depends on the microscopic param- 
eters of the problem. We have considered the most inter- 
esting case {Is <C II and d < a < lu ), when all possible 
regimes are developed. For other choices of parameters, 
some of the regimes may disappear (see, in particular, 
the footnote Q). 



IV. CONCLUSIONS 

In this paper, we have studied the Loschmidt echo (or, 
in a different terminology, the fldelity), which character- 
izes the sensitivity of a quantum system to an external 
perturbation, in a model with a weak quantum random 



potential. Using the diagrammatic approach, we have 
shown that at short times the fidelity decays exponen- 
tially with the rate 2/fs set by the golden rule, while its 
long-time asymptotics is of power-law type and is deter- 
mined by the diffusive nature of the dynamics on this 
time scale. For a sufficiently long-range disorder a time 
range emerges where the diagrammatics becomes inap- 
plicable due to merging of two ballistic diffusons into a 
more complicated "4-diffuson". To study the fidelity in 
this regime, we have applied a quasiclassical (path inte- 
gral) approach. This allowed us to express the disorder- 
averaged fidelity in terms of a solution of a quasiclassical 
evolution equation, see Eqs. (33) and (p4[). On time scales 
larger than this equation takes a form of the kinetic 
equation for the distribution function g{(j)^p-.) of phase- 
space separations between two classical paths (one of 
which is subject to the perturbation). Solving the kinetic 
equation, we find several additional regimes of behavior 
of the Loschmidt echo, separating the short-time golden- 
rule decay from the long-time diffusive asymptotics. In 
particular, there arises a "Lyapunov regime" , where the 
fidelity decays exponentially with a rate governed by the 
classical Lyapunov exponent. 

It is worth mentioning that our path-integral calcula- 
tion of the fidelity is closely connected to the analysis of 
quantum interference effects in a long-range disorder per- 
formed in [p^ In particular, after the Fourier trans- 
formation — !■ our evolution equation ( ^8| ) has the 
same form as the equation describing the "Hikami box" in 
]l9[ (after averaging over the smooth random potential). 
This is a remarkable agreement, since the methods used 
are essentially different: Aleiner and Larkin work in 
a given realization of a random potential (which assumes 
that the potential is classical, i.e. the condition oppo- 
site to Eq. (^) is fulfilled), while we consider the case of 
a diffractive scattering [Eq. (^] and perform all calcula- 
tions for disorder-averaged quantities. There is, however, 
an important difference between the equations obtained. 
Specifically, in our case the last term of Eq. ( ^8| ) (which 
is proportional to l/fn) is due to the difference 6U be- 
tween the Hamiltonians for the forward and the back- 
ward propagation. On the other hand, the authors of 
Ref. add "by hand" a term of exactly the same type 
(with a certain time Tq replacing our ftr) for a problem 
without any perturbation SU, arguing that it mimics a 
small-angle diffraction in the system. To our opinion, this 
justification is questionable (at least, for a system with 
a weak smooth disorder). Indeed, in this case all scat- 
tering processes determining the transport in the system 
are of diffractive type and are taken into account in our 
approach. There is thus no freedom to add an additional 
"diffractive" term to the kinetic equation. We thus be- 
lieve that the Hikami box is described by an equation 
without such term (i.e. analogous to our equations (|34|), 
( ^ ) in the absence of perturbation, SU,f^^ = 0) but 
with appropriate boundary conditions. While this will 
probably not affect the main results of Ref. (depend- 
ing only logarithmically on r^), such a more consistent 
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treatment of the quasiclassical Hikami box would 
be of conceptual importance for the theory of quantum 
interference effects in systems with large-scale inhomo- 
geneities. We leave this issue as an open problem for the 
future research. 
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